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THE ROLE OF TIHE-HI STORY EFFECTS IN THE FORMULATION 
OF THE AERODYNAMICS OF AIRCRAFT DYflAMICS 

Murray Tobak* and Lewis B. Schiff** 

Ames Research Center, NASA, Moffett Field, California 94035, U.S.A. 


SUMMARY 

The scope of any aerodynamic formulation proposing to embrace a range of possible maneuvers is shown 
to be determined principally by the extent to which the aerodynamic indicial response is allowed to 
depend on the past motion. Starting from the linearized formulation, in which the indicial response is 
independent of the past motion, two successively more comprehensive statements about the dependence on 
the past motion are assigned to the indicial response (1) dependence only on the recent past and (2) depen 
dence additionally on a characteristic feature of the distant past. The first enables the rational 
introduction of nonliriar effects and accommodates a description of the rate-dependent aerodynamic 
phenomena characteristic of airfoils in low-speed dynamic stall; the second permits a description of the 
double-valued aerodynamic behavior characteristic of certain kinds of aircraft stall. An aerodynamic 
formulation based on the second statement, automatically ernbracirg the first, may be sufficiently compre- 
hensive to include a large part of the aircraft's possible maneuvers. The results suggest a favorable 
conclusion regarding the role of dynamic stability experiments in flight dynamics studies. 

LIST OF SYMBOLS 

0^1 pitching-moment coefficient, P— - ^ ^ ^ 

G[o(c)] functional notation: value at C * t of a time- dependent .'unction which depends on all values 

taken by the argument function a(c) over the fi.me interval 0 < t < t 

I moment of inertia about the pitching a ns 

I reference length 

q dynamic pressure, y pV’ 

S reference area 

t time 

V magnitude of flight velocity vector 

0 angle of attack, Fig. 1 

p atmospheric density 

u frequency of harmonic oscillatory motion 

1. INTRODUCTION 

One of the difficult problems in aircraft flight dynamics is that of formulating an aerodynamic force 
and moment system with sufficient scope to cover the wide range of maneuvers typical of modern aircraft 
(Refs. 1,2). What is the nature of the problem? 

Consider the questions that arise in the prediction of a maneuver from a known Initial state. Let an 
essentially rigid aircraft with known inertial properties undergo an arbitrary motion. At a certain time 
to. allow a measurement of the aircraft's state (i.e., its linear and angular velocity components) and 
its aerodynamic response (i.e., the aerodynamic force and moment). Given this information at to. what is 
needed to predict the aircraft's motion over the succeeding Increment of time? The ability to carry the 
motion forward over the first Increment of time implies, of course, the ability to predict the entire 
subsequent motion. What is needed principally is a form for the incremental changes in the aerodynamic 
force and moment, that is, the indicial response, over the Increment of time. Assigning an adequate form 
constitutes the problem of formulation. The difficulty of the problem arises in assigning a form that 
applies not only to the motion under study, but to all of the other motions of which the aircraft is 
capable, and which might have occurred prior to tQ. This way of describing the difficulty allows one to 
appreciate the great virtue of a linearized version of the aerodynamic indicial response. Invoking 
linearity supposes that the aerodynamic indicial response is independent of anything that happened prior 
to the origin of the response. Thus, the calculation can be carried forward without any acknowledgment 
whatever of the motion prior to tp. Although there are flow regimes where use of the linearized formula- 
tion can be justified (e.g., attached flows with small perturbations), these regimes do not embrace the 
whole range of flows that a modem aircraft may experience. A fomulation applicable to the remaining 
regimes must be freed of the limitation Imposed by linearization. This means, of course, that the 
aerodynamic indicial response must be allowed to depend on the past motion. 

In a series of papers (cf. Ref. 3 for a connected account), the authors have tried to show how con- 
cepts from functional analysis could be used to construct a mathematical framework allowing a general 
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dependence of the aerodynamic Indiclal response on the past notion. Having a rigorous framework has 
enabled the Introduction of rational approximations which, in effect, limit the dependence on the past 
motion to some definite property, applicable to an appropriate class of flows. One can argue In favor of 
limited statements about the dependence on the past motion as follows: Since tne linearized formulation 

has found application to a certain class of flows, a formulation based on a limited statement, which 
Includes the linearized formulation as a special case, must find application to a wider class of flows. A 
sequence of successively more comprehensive statements, each emtuacing all of the preceding ones, must 
eventually reach a stage where the resulting formulation Is applicable to a sufficiently wide class of 
flows to accommodate a description of all of the .notions of Interest. It remains to discover whethp' this 
stage can be reached well short of having to account for the whole past motion In detail for any of the 
motions of Interest. 

Thus, the role assigned to time-history effects, that Is, the statement about the dependence on the 
past motion, constitutes a determining criterion by which the 'nerlts and shortcomings of any aerodynamic 
formulation may be judged. The purpose of this report Is to investigate how far the first few statements 
go toward fulfilling the goal of a sufflclentlv comprehensive statement. Starting from the linearized 
formulation. In which the Indiclal response is independent of the past motion, two successl'^ly more com- 
prehensive statements about the dependence on the past motion are assigned to the Indiclal response 
(1) dependence only on the recent past and (2) dependence additionally on a characteristic feature of the 
distant past. The successive statements allow the effects of successively larger bodies of aerodynamic 
phenomena to be acknowledged within the scopes o^ the resulting formulations. The first enables the 
rational Introduction of nonlinear effects and acconmodates a description of the rate-dependent aerodynamic 
phenomena characteristic of airfoils In low-speed dynamic stall; the second permits a description of the 
double-valued aerodynamic behavior characteristic of certa’n kinds of aircraft stall. It Is suggested that 
an aerodynamic formulation based on the second statement, automatically including the first, may be of 
sufficient scope to embrace a large part of the aircraft's possible maneuvers. Implications of the results 
with regard to dynamic stability experiments are discussed. 

2. DEFINITION OF MANEUVER 

To focus directly on the question of time-history effects, It is advisable to avoid the complications 
Introduced by coordinate systems and motions with multiple degrees of freedom. In all of the study to 
follow, the aircraft's maneuver Is restricted to be planar with only a single degree of freedom. Extension 
to more general motions will be straightforward, paialleling that described in Ref. 3. 

Let the aircraft be In level steady flight prior to time zero. At time zero let It begin an 
arbitrary pitching maneuver du'^ing wnich the center qravUy continues to follow a rectilinear path at 
constant velocity V. Hence flight-path properties such as dynamic pressure, Mach number, and Reynolds 
number remain constant throughout the maneuver. The pitching maneuver is defined by the angle of attack a 
(Fig. 1), the angle between the aircraft's longitudinal axis and the velocity vector. The motion, of 
course, may be specified to reproduce that of a wind-tunnel mjdel in an osclllatlons-in-pitch experiment. 
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Fig. 1. S1ngle-degree-of-f reedom pitching maneuver. 


Thus, focusing on this motion will facilitate a later discussion of the implications of the results with 
regard to wind-tunnel experiments. 

3. FORMATION OF INDICIAL PITCHING-MOMENT RESPONSE AND INTEGRAL FORMS 



Since it will be necessary In later 
sections to consider the Influence of random 
fluctuations, the formation of the Indiclal 
pitching-moment response will be described In 
a way that acknowledges their presence. Two 
motions have to be considered (cf. Fig. 2). 
First, beginning at t * 0, the aircraft Is 
made to undergo the motion under study o(c). 
At a certain time t the motion Is con- 
strained so that the value of o at time t, 
that Is, oh), remains constant thereafter. 

The pitching moment corresponding to this 
motion is measured at a ^Ime t. Now If o(t) 


Fig. 2. Formation of Indiclal response. 
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is Sufficiently large so that, fo, example. f1cv< separation oirurs in ttie couric of a maneuver, then ni a 
result of the ensuing fluctuations in the flo*y, an> single measurement of the pitching moment at time t 
will Include a random component. This circumstance calls for repeating the maneuver and the correspond! m, 
measurement at time t 'iiany times and taking the arithmetic mean of the measurements. If the fluctuating 
part of the response is truly random, its contribution to the measurement at time t should cancel in 
the mean, and the resulting mean value should be representative of the deterministic part of the response. 
It will be assumed that this is true for any time t, and that as a result, a deterministic part cf the 
response will exist that 1$ continuous for all c i» the interval 0 < t < t. Second, the aircraft is 
made to execute precisely the sane motion, beginning at c ° 0 and constrained in the same way at c * . 

except that at the latter time, i is given an incremental step to over its value at t - t. Hentt, for 
all time subsequent to t, is equal to a(r) ♦ ac. The pitching monent corresponding to tnis irotiou 
is again measured at time t. Just as before, the second maneuver and the corresponding measurement at 
time t must be repeated many times, and the arithmetic mean of the measurements taken to be the deter 
ministic part of the response. The difference between mean values for the two measurements. aC^{t), is 
divided by the incremental step ao. The limit of this ratio (if it exists) as the magnitude of the step 
approaches zero Is called the indicial pitchinq-moment response at time t per unit step change in at 
time T. Since the two maneuvers prior to t = t are identical (in the mean), the ratio must be identi- 
cally zero for 0 < c < t. At c = t a discontinuity in the ratio is permissible, reflecting the 
discontinuous change in o. For all c > t the ratio is assumed to be continuous. With the understandiii 
that the pitchinq-moment response to each maneuver and at each time t is the result of an ensemble 
average of measurements, the Indicial pitching-moment response is defined as: 


1 im 

ao-*-o 




Cm [o(5);t,i] 
c 


( 1 ) 


As the functional notation indicates, the indicial response is allowed to depend in an unspecified way on 
the entire history of the irotion o(;). 


When the asswptions leading to the definition of a deterministic indicial response can be said to 
hold within each increment of the stepwise representation of an arbitrary motion o(t), the pitchinq- 
moment response Cm(t) to the motion o(t) follows from a summation of incremental responses over the 
time interval 0 to t: 




= Cm(0) 


•/: 


Cm [o(0;t.T] 


( 2 ) 


This is the general integral form for Cq,(t) corresponding to an arbitrary motion o(t). The form is 
essentially exact, but its further use without approximation is exceedingly difficult. The nature of the 
difficulty becomes clear if one writes the equation of motion for the si nqle-degree-of- freedom pitching 
motion a(t), and asks for a solution of a(t) for specified initial conditions. The equation of motion 
is: 


o(t) - (3f^) |c„(0l +jr^Cm^(^(t);t,T] ^ di} 


<3' 


Since the indicial response within the inteqral is a functional, dependent in general on the whole past 
motion a{t). it is unknown when o is unknown. Thus, both the indicial response and the motion must be 
found simultaneously, an awesome prospect. Cases can be envisioned (e.q., massively ablating reentry 
vehicles) where the mutual dependence between the motion and the indicial response cannot be uncoupled. 
However, for rigid aircraft of fixed shape, the past success of the linearized formulation lends credence 
to the belief that the interdependence may be at least partially uncoupled. Now if the indicial response 
were somehow known. Eg. (3) would become an inteqro-differential equation of the Volterra type for which 
solutions can be found by known, albeit numerical, techniques. The indicial response could be considered 
known, at least in principle, if its dependence ,m the past motion were specified In a way that allowed ii 
to be an identifiable member of a collection of indicial responses, all of which had been obtained btfore- 
hand from, for example, a suitable series of experiments. Thus, Eg. (3) can be made tractable by assignin 
to the indicial response .’.ppropriaie statements about how it depends on the past motion which allow it to 
be determined in advance. Let it be noted, however, that since every statement assigned to the indicial 
response will multiply the number of responses in the collection required to be known in advance, it 
becomes imperative to make the least number of statements possible. In succeeding sections, two such 
statements will be introduced, which it is hoped, may suffice to cover the cases of interest. 

Finally, it is recognized that counting on the availability of a collection of indicial responses 
may be somewhat unrealistic in view of the difficulty of experimentally determining an indicial response. 
In a later section, however. It will be shown how the results from suitably designed oscillations-in-pitci. 
experiments, which. It may be presumed with greater reason, are technically feasible, may be used in plai.e 
of the integral In Eg. (3). Thus, for the availability in principle of a collection of Indicial response' 
may be read the availability of an equivalent collection of results from oscil lations-in-pitch experiments 


4. DEPENDENCE OF THE INDICIAL RESPONSE ON THE RECENT PAST 

First, £o. (1) for the indicial pitchinn-moment response will be put in an equivalent form that wi'1 
suggest a first statement about its dependence on the past motion. If ci(r) can be considered to he 
analytic in a neighborhood of f, • t (corresponding to the most recent nast for an indicial response with 
origin at t • t), in principle, its history can be reconstructed from a knowledge of all of the coeffi- 
cients of its Taylor series expansion about t * '• Then, since n(f.) is equally represented by the 
coefficients of its expansion, the functional, with its dependence on o(c), can be replaced without 
approximation by a function with a dependence on all of the coefficients of the expansion of •(;) about 
t • t; that is. 
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(t - T;o(T).J(T),'j(T), . . .) (4) 

o "o 

The additional replacement of a dependence on el apsed time t - t rather than on t and t separately Is 
Justified within the soeclficatlon already Invoiced of constant fllqht-path properties. 

Now It Is argued that a class of flows exists for which disturbances originating at times far removed 
from the vicinity of c • t will have died out before they are able to Influence events In the vicinity 
of c • T. In such cases. It can be assumed that the Indiclal response will have "forgotten" long-past 
events, and so will depend only on events occurring in the most recent past. Therefore, to the extent that 
the indiclal response can be influenced by the past motion, the form of the past motion just prior to the 
origin of the step might just as well have existed for all earlier times. Accordingly, only the first few 
coefficients of the expansion of o(c) need be retained to characterize correctly the most recent past, 
which Is all the indiclal response is assumed to remember. Retaining the first two coefficients of a(c), 
for example, implies matching the true past history of o in magnitude and slope at the origin of the 
step, thereby approximating o(t) by a linear function of time 3 (c) = o(t) - 4(t)(t - C). With an 
approximation of this order In force In Eg. (1), the integral form replacing Eg. (2) and the right-hand 
side of Eg. (3) becomes 

C^(t) = C^(0) C„,^(t - t;o{t).5(t 1)^ dr (5) 

The steady-state value of the iridiclal response can be put In evidence with the additional (consistent) 
assumption that events In the recent past, that Is for c < t, will again be far removed and so, forgotten, 
so far as the Indiclal response Is concerned when t - t - This means that the steady-state value of 
the Indiclal response will depend only on local conditions, that is, on the constant value of o(t). The 
latter behavior Is put in evidence by the substitution 

(t - t;o(t),5(t))= {“';o(t)) - r(t - r;c(T),4{T)) (6) 

0 3 

where Cm,(‘*;o(T )) Is the steady-state value of the Indiclal response. Notice that It must be a single- 
val ued function of o(t). The function F(t - ’■:p(t),J{t)) Is called the deficiency function; it approaches 
zero as t - T - and. In practice, will be essentially zero for all elapsed time t - t larger than a 
relatively small value t*. When En. (6) Is substituted in Eg. (5), the steady-state term multiplied by 
(do/drldr forms a perfect differential which can be Integrated. The resulting formulation for Cmlt) 
becomes 

C^(t) = C|^(";r(t)) - J P(t - T;o(T),j(t)) dt (7) 

Here, C^(»;o(t)) Is the pitching-moment coefficient that would be measured in a steady flow with 3 fixed 

at the instantaneous value o(t). Again, note that ^6 a single-v a lued function of o 

according to this formulation. 

Eguation (7) actually Includes three Increasingly comprehensive formulations, each of which may be 
applicable In appropriate circumstances. The simplest, of course, is the linear formulation, for which 
the indiclal response Is said to be Independent of both o(x) and 4 {t). The resulting simplification 
Is reflected in the eguation of motion. Eg. (3), wMch becomes 

ij(t) = -J F(t ' t) ^ (8) 

The eguation Is linear and the Integral term Is of the convolution type, which enables an iimedlate 
solution for o by the aid of Laplace transforms. A considerable additional virtue Is that the collec- 
tion of Indiclal responses regulred to be known In advance consists of one member. 

The second formulation, anpllcable In particular to slowly varying motions. Is obtained by omitting 
the dependence on 4(x) from the Indiclal response. Omitting this dependence In Eg. (6) means that, so 
far as the indiclal response Is concerned, the motion prior to the origin of a step Is being aoproxlmated 
by the time-invariant motion a(r,) . o{x). The equation of motion. Eg. (3), becomes 

a nonlinear Volterra Integro-differentlal eguation, solvable by numerical techniques. Here, the collection 
of Indiclal responses required to be known In advance must consist of members corresponding to a range of 
values of o{r). This formulation also lends Itself to reduction to a form correct to the first order In 
frequency, resulting In a nonlinear generalization of the classical stability derivative formulation. The 
form has been studied at length In the authors' previous work (Ref. 3). Approximation at the level of the 
second formulation thus enables the rational Introduction of nonlinear effects. 

Finally, the third, most comprehensive formulation Is that represented In full by Eg. (7). This form 
Is of sufficient scope to allow the treatment of motions Involving hysteresis effects caused by rate- 
dependent aerodynamic phenomena. It Is believed to be applicable, for example, to the complex set of 
aerodynamic phenomena characteristic of airfoils In low-speed dynamic stall (Refs. 4-6). Retaining a 
dependence on n(x) allows assigning different Indiclal responses to a step at a single value of o(x), 
depending on the magnitude and sign of o(x). It Is possible, for example, to distinguish between Indiclal 
responses where o was Increasing or decreasing prior to the step. The motion that would be required to 
oht.iln the Indldal response experimentally with positive 6(x) Is essentially the same as the maneuver 
that has been used fo study the overshoot In lift that occurs following a rapid pitch-up (Ref. 7). It would 
alto he netPssarv tn carry out experiments Involving pitch-down maneuvers to allow for the possibility 
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that the Indicia! response for a given o(t) and a negative value of o(t) will be different from that 
for a positive 6(t) of the same magnitude. Here, the collection of Indlclal responses required to be 
known In advance must embrace a range of values of <i(t) and a range of both positive and negative values 
of 5(t). 

While Eq. (7) goes some distance toward fulfilling the goal of a sufficiently comprehensive 
formulation, it Is still Incapable of accounting for the existence of multivalued aerodynamic phenomena 
that need not depend on the pitching rate. This Is evident, In particular. In the representation of the 
steady-state aerodynamic pitching moment Cn,(«*;o(t)), and Its derivative Cmg(«;!>(t) ) , which must be 
single-valued functions of o. Admitting the possibility of multivalued aerodynamic responses, not 
necessarily dependent on the pitching rate, will require an acknowledgment of the Influence of the distant 
past on the Indlclal response. 


5. DEPENDENCE OF THE INDICIAL RESPONSE ON THE DISTANT PAST 

In the preceding treatment, on the assumption that events In the distant past should be Incapable 
of Influencing the Indlclal response, the Indlclal response functional Cmg[o(c);t,T] was replaced by a 
function CiH(j(t - t;o(t),o(t)) which depends only on the magnitude and slope of the past motion c(c) 
at the origin of the step. This replacement can be viewed either as an approximation of the actual motion 
o(c) by a linearly varying motion o(c) » o(i) - *(t)(t - c) fo'" past time, or, as a substitution 
applicable only In the vicinity of f. » t with the imnllcit understanding that the distant-past motion, 
being Immaterial to the Indlclal response, can be assigned at will. In what follows, the latter Inter- 
pretation will be the desired one to the extent possible. However, the presence of fluctuations appears 
to represent a condition where events occurring In the distant past (e.g., the initiation of fluctuations 
due to flow separation) can affect the present evolution of the Indlclal response. For suppose that 
certain distant-oast motions, but not others, can Initiate fluctuations that persist up to the measuring 
time, and that the Indlclal response measured 1n the presence of fluctuations can differ from the measure- 
ment In the absence of fluctuations. Then this must be taken Into account somehow, and In a way that 
does not require assigning the ac ual distant-past motion to every indlclal response. 


5.1 Reformulation of Pitching-Moment Response 


It Is convenient to consider first a typical motion qualifying as the first of the two motions 
required to form the Indlclal response (cf. Fig. 3). The motion for values of c, < t Is given by o(t); 
It Is held constant at o{t) for all c > t. 

The time c ' t Is the time at which the 
pitching moment Is measured. The whole 
time period prior to c • t will be 
called the past relative to any value 
of c ” t > T. Let the past be divided 
into two parts so that there 1s an 
Interval T just prior to E * t. 

The Interval of duration T will be 
called the recent past and the 
remaining interval the distant past . 

Let T be chosen sufficiently large 
so that under normal circumstances, 
the particular form of the motion 
o(E) in the distant past Is Inmate- 
rial to the pitchinq-moment response 
for values of ^ * t > t. This is 
the condition that has been used in 
the preceding section. 



Fig. 3. Recent past and distant past. 


How can motion In the distant past Influence the measurement for the pitching moment at time t? 

It Is clear that It will do so if it Is capable of changing the nature of the flow in the recent past, 
to which the measurement certainly sensitive. This is possible if motion in the distant past has 
initiated fluctuations that persist into the recent past. The following Is an explicit argument. 

Suppose that the pitching moment is measured at time t for the motion shown in Fig. 3. (Recall that 
the measurement Is dCtually the ensemble average of repeated measurements at t corresponding to 
repetitions of the same motion.) With everything else remaining the same, make an Infinitesimal change 
in the distant-past motion as shown by the dotted curve In Fig. 3. After many repetitions of this motion 
and the corresponding measurement for the pitching moment at time t, compare the ensemble-averaged 
result with that of the original motion. Suppose first that there is no difference between the two 
measurements. This, of course. Is the expected result according to the preceding formulation. However, 
there are two conditions under which It should hold: (1) neither distant-past motion has Initiated 
fluctuations that persist Into the recent past or (2) both distant-past motions have Initiated fluctuations 
that attain similar statistical properties over the identical recent past, thereby altering the ensemble- 
averaged flow over the recent-past motion In the same way. Both conditloni can be Incorporated within 
a single assertion that In neither case did only one motion Initiate fluctuations persisting Into the 
recent past. Now suppose that the two measurements differ at time t. Then the Infinitesimal change 
In the distant-past motion must have exceeded a cri tical condition for the Initiation of fluctuations, 
and the persistence of these fluctuations into tHF recent past must have changed the form of the flow In 
the recent past In ensemble average. The argument of exceeding a critical condition Is equally valid for 
cases In which an Infinitesimal change In the distant-past motion stops previously existing fluctuations 
from persisting Into the recent past. Thus, If the distant-past motion Is able to Influence the measure- 
ment at t, the cause Is attributable to the existence of a critical condition for the Initiation or cessa- 
tion of persistent fluctuations. 
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Finally, consider a second Infinitesimal change In the distant-past motion that also exceeds the 
critical condition. Consistent with condition (?) already noted, It Is argued that the pitching-moment 
response to this motion at time t will be the same as the response to the previous motion which had 
also exceeded the critical condition. The argument Is that relative to values of c ■ t > t, there Is 
sufficient time over the duration T of the recent past for fluctuations originating In the distant past 
to attain statistical properties that no longer depend on their origin , but rather depend only on their 
experience over the recent past. Since the recent-past motion Is tne same for both distant-past motions, 
the statistical properties of their respective fluctuations should have become the same by the time they 
have reached the vicinity of c ■ t. With Identical statistical properties, ensemble averaging of 
repetitions of the pitching-moment responses to the two motions should then yield the same result at ( • t. 
Now the same argument should hold for an^ distant-past motion that has Initiated fluctuations persisting 
Into the recent past. According to this argument then, the pitching-moment response at time t for a 
given past motion can be duplicated by the pitching-moment response for a motion whose form In the distant 
past Is assigned at will, so long as It Is known that both have Initiated fluctuations persisting Into a 
same recent past. Likewise, a distant-past motion whose fluctuations do not persist Into the recent oast 
can be replaced, according to the argument, by one assigned at will, so long as fluctuations originating 
In the latter motion also do not persist Into the recent past. The argument translates into a mathematical 
statement that the pltchlng-momerit response at c « t > t must be a unique functional of the recent-past 
motion and additionally, must depend on a parameter that designates by say, one of two numbers, whether 
fluctuations originating In an otherwise arbitrary distant-past motion do or do not persist Ifito the recent 
oast. Accordingly, the pitching-moment response Is written In the functional form 

-ni(t) • C„,[o(t);t,T,»(T)] (10) 

where It will be understood that the range of c over which there Is a functional dependence on o(c) Is 
now restricted to cover the recent past only, that Is, t - T < c < t, and where X designates the type 
of distant-past motion, otherwise arbitrary, by one of two numbers; that Is, 


x(t) • 0: fluctuations originating In distant past do not persist Into recent past; 

x(t) ■ 1: fluctuations originating In distant past persist Into recent past. 


( 11 ) 


Notice In Eqs (10) and (11) that although X can take only one of two values, nevertheless It Is 
denoted a function of t. This Is because the event characterized by x must originate In the distant 
past, and by definition, what Is called the distant past depends on t. That Is, an event Is said to occur 
In the distant past If the time t at which It occurs satisfies t - t > T. 


5.? Reformulation of Indiclal Response and Integral Form 


Now It Is necessary to consider the role of fluctuations In the formation of the Indiclal response. 
Consider again the motion Illustrated In Fig. 3, for which the pitching-moment response Is of the form 
Eg. (10). Just as before, with everyth!, g else remaining the same, make a small change ao In the con- 
stant value o(t). Repetitions of this motion and the corresponding measurement at time t yield an 
ensemble-average value 'or the pitching moment at time t. The procedure Is repeated for successively 
smaller values of Ao, as many times as necessary, to carry out the limiting process Indicated formally 
In Fq. (12): 


11m 

AO-*0 


Ac 


11 

Ao-»0 


Im 

-»o t 


[oj(c);t,T,v(T)] - 




■,x(t)] 


An 


* Cnt [o(t);t,T,XlT)] (12) 

a 

where 

o,(t) * o(t) ; 0 «. C < T 

* <j(t) ; C > I 

02(e) « o(e) ; 0 < e < T 

*o(t)*Ao; e^T 


The result (If It exists) Is the definition of the Indiclal response. Suppose first that the limit 
exists, and moreover, that the pitching-moment response to the second motion Is very similar In character 
to that of the first. Again, there are two conditions; (1) neither motion oi(t) nor 01 (c) contains 
flow fluctuations or (2) both motions 0 |(t) and 02 (c) contain flow fluctuations of similar statistical 
character. Both conditions are Incorporated within a single assertion that In neither case did the 
Infinitesimal change Ao either Initiate fluctuations or stop prior fluctuations. Now suppose that the 
ensemble-averaged pitching-moment response to the second motion differs from that of the first In such a 
marked way that the limiting process (ACfl,(t)/Ao) almost falls to converge. In this event. It will be 

said that the Infinitesimal change Ao has either Initiated fluctuations or stopped prior fluctuations. 

In other words. If the Infinitesimal change Ao required to fonti the Indiclal response Itself either 
Initiates fluctuations or stops prior fluctuations, the evidence of this will be a marked, almost 
discontinuous, change in the Indiclal response. This abrupt change will be evidenced. In particular. In 
the steady-state va’ue of the Indiclal response. While the cause of an abrupt change In the Indiclal 
response most o^ten can be attributed to the Initiation or termination of fluctuations, let It be noted 
that nothing In the analysis prevents associating the abrupt change with other, perhaps nonfluctuating. 
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phenomena as well; for example, an abrupt shift In the pattern of flow separation or an abrupt shift In 
the location of a shock wave (Ref. 8). All that the analysis requires In principle Is the existence of 
two distinct reqlmes of flow separated by critical conditions. Thus, the phrase "Initiation (or termina- 
tion) of fluctuations" may stand for any flow phenomenon leading to an abrupt change In the Indiclal 
response. 

To complete the reformulation. It Is necessary to consider the suimatlon process leading to an 
Integral form for the pitching-moment response to an arbitrary motion. Although, as noted, the possibility 
Is allowed of a near-dlscontlnul ty with respect to a In the Indiclal response, this possibility will exist 
only at certain Isolated values of o where the exceeding of a critical condition by the Infinitesimal 
change so either Initiates or stops fluctuations. Since the nearly singular behavior of the Indirlal 
re'.ponse thus will be confined to discrete events In the history of o, these events will not Invalidate 
the general applicability of the sunmatlnn procedure. Then, as before, the result of sumiatlon yields for 
Cm(t): ^ 

‘ * I ['»(«);t,T,x(T)] ^ dT (13) 

•o ° 

The form differs from the one presented earlier (Eq. (2)) In two respects: (1) the functional dependence 

jf the indiclal response on o(t) extends only over the recent past; (2) the Indiclal response depends 
additionally on the parameter x, designating the type of distant-past motion that Is to be attached 
to the recent-past motion. 

Finally, the same argument that was used before can be Invoked to replace the Indiclal response 
tur.rtlonal in Eg. (13) by a function dependent on a limited number o^ oarameters, rather than on all values 
cf T(r.'. over the Interval T of the recent past. If it Is assianed again that the recent-past motion Is 
ddri'uatelv represented by the first two terms of Its Taylor series expansion about t * x, then Eq. (13) 
becomes 

^ dx (14) 

Again, It can be argued that as the Indiclal response approaches Its steady-state value with Increasing 
va’ues of t - t, It must become Independent of any particular recent-past motion (characterized by ?(t) 

'.n Eq. (14)), since the statistical properties of disturbances originating In the recent past will have 
becore Independent of their origin as t - t -* ■. The dependence on x remains, however, and this means 
th<>x ie steady-state value of the indiclal response may now be a dou ble-valued function of o(x), 
correspjnding to the two possible values ot x. The substitution 

Cm (t - t;o(t),o(t),x(t)) * Cm (";o(x) .X(x)) - F(t - t;o(t),o(t),X(t)) (15) 

0 ^ 

puts the stea ly-state value of the Indiclal response In evidence. If Eq. (15) Is substituted In Eq. (14) 
and the interval 0 < x < t Is divided Into segments, each of which contains only a single value of x, 
then the steady-state term can be Integrated over each of the segments. The Intermediate terms always 
cancel, however, so that the formulation for Cm(t) takes the form: 

C^(t) • Cm(";o(t).x(t))- F(t - x;o(x),o(x),x(x)) ^ dx (16) 

As before, Cm(»;9(t),' (t)) Is the pitching -moment coefficient that would be measured In a steady flow with 
0 Tlxed at the Instantaneous value o(t). Like Its derivative. It may now be a double-valued function of 
f corresoondlng to the two possible values of >. 


5.3 Decision Logic for Choice of x 

It remains to determine a logic for assigning the appropriate value of X to an Indiclal response. 
Recall that the influence of the distant-past motion on an Indiclal response at current time, characterized 
by X, was associated with the exceeding of a critical condition In the distant past which either Initiated 
or terminated persistent fluctuations. It was noted also that at the time of exceeding a critical condition 
the initiation or cessation of fluctuations would be evidenced by a marked change In the behavior of the 
Indiclal response, and In particular the steady-state value of the Indiclal response, accompanying an 
Infinitesimal change In p. This suggests that a suitable experiment for determining where the onset and 
cessation of fluctuations occur Is simply the experiment that would be required In any case to determine 
the steady-state pitching moment as a function of angle of attack. A suitable program for changing the 
angle of attack Is Illustrated In Fig. 4, where, at each level, sufficient time is allowed before measuring 
the pitching moment for the pitching-moment response 
to reach a steady state (In the mean). Advancing 
through a series of Increasing angles of attack, 
and then similarly, through a series of decreasing 
angles of attack, should allow determining 
whether double-val ued behavior of the steady- 
state pitching moment is possible. Suppose 
that the result of measurements for the steady- 
state pitching-moment coefficient resentles 
that shown on Fig. 5(a). The curve has two 
distinct branches, reflecting the nonfluctu- 
ating and fluctuating regimes of flow. A 
region of overlap, ag < o * 05, exists In which 
the pitching-moment coefficient Is double- 
valued. As Illustrated schematically In 
Fig. 5(b), passage from one regime to the 



Fig. 4. 


Anqle-of-attack program for measuring 
steady-state pitching moment. 
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(a) Pltchlnq-moment coefficient. {l4 dcglnes of flows. 

Fig. 5. Steady-state pitching-moment coefficient with region of double-valued dependence on angle 

of attack. 


other Is barred (Indicated by the double line) except by traversing the critical points u = us, o ■ cp 
in the directions indicated by the arrows. Then, having the values of and op should suffice to 
determine a logic for the choice of A. 

Assume that a step-by-step calculation Is being made of a maneuver, and that the calculation has 

the calculation one n;ore step, U Is necessary 
to assign the appropriate indicial response to 
the point . » I. Can it now le done? The form 
of the ino.-ial r'^'sponse, Fg (15), indi otes tut 
parameters that have to Le known: o(t), o(t), and 

a(t). Since o(c) is known for t < t, the values 
of o(t), o(t) can be specified. It remains to 
determine whether a » n or 1 to complete the 
specification. The followinn three question-, are 
asked: 0) Is there at least one (n with 

T - Co > T such that o(Cc) * (2) Is there 

at least one Ci with i - Ci > T such that 
o(Ci) * op? (3) If the answers to (1) and (?) arc 
yes, is m1n(T - C]) > m1n(t - Cq) 2 fes or no 
answers to the three guestlons determine the value 
of >. Results are given 1n Table 1: 


liRLF 1. OFCISinN LOf.IC FOR A 
1 ? 3 ~A 

No 0 

Yes No -- 1 

Yes Yes No 0 

Ye s Yes Yes 1 


The formulation Eg. (16) Is considered to be the principal resuU of this study. Given the informa- 
tion reguired to determine A, its scope may be sufficiently wide to embrace motions Involving both 
nonlinear aerodynamic responses and the double-valued aerodynamic behavior evident In certain kinds of 
aircraft stall (Refs. 8-10). However, its usefulness In practice appears to hinge on the availability 
of a collection of indicial responses, forming the kernel of the integral term in Eg. (16). As noted 
earlier, the eventual availability of such a collection is very unlikely in view of the great difficulty 
of experimentally determining an indicial response. The purpose of this section will be to Show that the 
usefulness of the formulation In fact Is not contingent on the availability of Indicial responses. It 
will be shown that the integral term in Eg. (16) is replaceable by results from a technically more feasible 
experiment: the osci 1 latlons-ln-pi tch experiment. 

For an arbitrary motion o(t), the contribution to Cn,(t) of the Integral term in Eg. (16) may be 
approximated as a finite sum of responses to discrete steps Ao(t). The form of the sumation, that is, 
the integrand, is illustrated in Fig. ?. At each step, the deficiency function F(t - a ;o ( t )A (t ) ,A (t 1 ) 
dies out to zero as t - r ■• «, and will be essentially zero for all t - : larger than a relatively 

small value t^. This is shown schematically In Fig. 7. It Is clear that in the surjriatlon of responses 

at time t, only the responses In the interval t - tj * t < t yield measurable contributions at t. 

Thus, the form of the motion o(r) outside the interval t - tj < t < t is immaterial to the sunmatlon, 

except insofar as it determines the value of A. On the assunption that the motion outside the interval 
is such that it ensures the correct value of a, the motion within the interval t - t. < i « t can be 
approximated without serious error by any convenient substitute motion. If, for example, the substitute 
notion Is harmonic, the resulting suiwiation of responses at time t will be that corresponding to an 
cguivalent harmonic oscillatory motion. 


advanced to the point £ « t (cf. Fig. 6). To continue 



Fig. 6. Step-by-step calculation of a maneuver. 
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Fig. 7. Sumratlon of Indiclal responses. 


An arbitrary harmonic motion 0 ^( 1 ) about a constant mean provides three arbitrary constants 
(assuminn that the frequency u Is chosen to match a characteristic frequency of the actual motion) 
which may be chosen to match three prooertles of the actual motion. It Is especially Important to match 
the actual motion In the Inmedlate vicinity of r • t where the contributions to the sunmatlcn are 
largest. If It Is chosen to match the actual va’'.ies of e{t), j(t), and cj(t) by a harmonic motion 
on('). then the form of the harmonic motion Is qlven by 

o^It) * [ait) ♦ sin «(t . ^ cos .(t - t) (17) 

su n, t 

6.(t) • i(t) cos .(t - t) - sin «(t - t) (18) 

n ji) 


SuLS’ltutlng Eqs. (17) and (18) for o(t), o(t) In the Integral term In Eq. (16) yields *or the Integral 
term (letting the lower limit be t - tg for consistency): 


J( 


t) • -o(t)Jj_j^ F(t - T;e^(t),o^(T),i(T))cos .-(t - ■t)di 


d(t 


t-t. 


F(t - T;?>,(T),4^(T),A(t))s1n u>(t - t)di 


(19) 


or, with a change of variable. 

J(t) • -oft) J ^ F(u;e^(t - u),o^(t - u),»(t - u))cos uiU du 


Hilf" f(.-„ 


ft - u),5^(t - u),i(t - u))s1n jU du 


( 20 ) 


This Is the contribution to Cm(t) from the Integral term, and It Is the same contribution that would be 
obtained from an osc111at1ons-(n-o1tch experiment for an oscillation constructed according to Eg. (17). 

This means matching o(t), o(t), and 3(t) requires that an equivalent harmonic motion have a mean value equal 
to [o(t) ♦ o(t)/u’] and an amplitude equal to ([5(t)/*]^ ♦ [3(t)/w‘]^)’/ The contribution from the 
term multiplied by ci(t) In Eq. (20) Is actually of second order .n w and probably negligible for the 
very low reduced frequencies typical of most aircraft motions. For rapid maneuvers, however, (e.g., a 
rapid pitch-up maneuver) It may be necessary to retain the term, since It accounts for the very large 
deviation from the steady-state aerodynamic contribution (I.e., Cn,(-;o(t),»(t)) In Eq. (16) at 0 • 

(where c E 0) which has been observed In wind-tunnel experiments with oscillating airfoils In the low- 
speed dynamic stall rejime (Refs. 4,5). 

The result repres'.-nted by Eq. (20) mathematical ly expresses the major theme of this study. However, 
couched In more general terms. It leads to an almost self-evident conclusion that should hold In the 
general case' the instantaneous aerodynamic force and moment corresponding to an arbitrary motion can be 
duplicated with the Instantaneous force and moment corresponding to an assigned motion, so long as both 
motions are essentially the same In tne recent past relative to the instant, and are In the same flow 
regime determined by the otherwise Immaterial distant-past motion. 


7. CONCLIAMNG RE14ARKS 

The scope of any aerodynamic formulation proooslnq to embrace a range of possible aircraft maneuvers 
has been shown to be determined principally bv the extent to which the aerodynamic Indiclal response Is 
allowed to depend on the past motion. Allowing the Indiclal response to depend onlv on motion In the 
recent past resulted <n an aerodynamic formulation enabling the rational Introduction of nonlinear effects 
and a descriotlon of the rate-dependent aerodynamic phenomena characteristic of airfoils In low-speed 
dynamic stall. Allowing the Indiclal response to de:>end additionally on a characteristic feature of 
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motion In the distent past, that Is, the Initiation or tennlnatlon of persistent fluctuations, resulted 
In a more comprehensive formulation permlttlnq a description of the double-valued aerodynamic behavior 
characteristic of certain kinds of aircraft stall. The scope the latter formulation should be 
sufficiently wide to Include any pitching maneuver having no more than two distinct regimes of flow, 
separated by critical conditions. Straightforward e<tens1ons of the formalism already developed should 
yield formulations permuting a description of any number of flow regimes and embracing motions with 
multiple degrees of freedon.. A general conclusion that can be drawn from this study, favorable regarding 
the role of dynamic stability experiments. Is the following: the Instantaneous aerodynamic force and 
moment coi responding to an arbitrary motion can be duplicated with the Instantaneous forte and moment 
corresponding to an assigned motion, so long as both motions are essentially the same In the recent past 
relative to the Instant, and are In the same flow reginie determined by the otherwise Inmaterlal distant- 
past motion. 
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